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SUMMARY 


This paper will summarize recent work that has been performed under 
a NASA Langley Contract NAS1-17864, "Additional Development of the XTRAN3S 
Computer Program." Work on this code, for steady and unsteady aerodynamic 
and aeroelastic analysis in the transonic flow regime, has concentrated in 
four areas: 

1) Maintenance of the XTRAN3S code, including correction of errors, 
enhancement of operational capability, and installation on the 
Cray X-MP system at Wright-Patterson Air Force Base; 

2) Extension of the vectorization concepts in XTRAN3S to include 
additional areas of the code for improved execution speed; 

3) Modification of the XTRAN3S algorithm for improved numerical 
stability for swept, tapered wing cases and improved computa- 
tional efficiency; 

4) Extension of the wing-only version of XTRAN3S to include pylon 
and nacelle or external store capability. 


• MISCELLANEOUS PROGRAM CORRECTIONS & ENHANCEMENTS 

• IMPROVED VECTORIZATION OF XTRAN3S 

• ALGORITHM IMPROVEMENT 

• DEVELOPMENT OF WING/PYLON/STORE CAPABILITY 
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MISCELLANEOUS PROGRAM CORRECTIONS 


This figure lists the various UPDATE correction sets for correction of errors 
in the XTRAN3S code. These correction sets are to be applied to BCS 
Version 1.6 (equivalent to Langley Version 1.5) of the XTRAN3S code. 


UPDATE NAME FUNCTION 

EXT LOFT CORRECTS LOFTING PROCEDURE FOR EXTRAPOLATION 

OF AIRFOIL DATA 


VISCORR 

FIXPROMAT 

MISCFIX 

ENCODE 

VECTCORR 

PARMGEO 


CORRECTS ERROR IN WEDGE PARAMETER DEFINITION 

CORRECTS DATA INPUT ERROR FOR INERTIA MATRIX 

CORRECTS GEOMETRY DEFINITION ERROR, MODIFIES 
CPRESS AND CFOMO TO ACCOUNT FOR MODIFIED GRID, 
AND VARIABLE DIMENSION ERROR IN ZSWEEP 

CORRECTS LABELING ERROR IN TIME HISTORY OUTPUT 

CORRECTS INITIALIZATION ERROR IN VECTORIZATION 

PROVIDES VARIABLE DIMENSION FOR NUMBER OF 
POINTS ON AIRFOIL 


MISCELLANEOUS PROGRAM ENHANCEMENTS 

This and the following figure list the UPDATE correction sets for 
improvement of XTRAN3S operational features. This group of correction 
sets, plus the previous list of error correction modifications, constitute 
upgrade of the XTRAN3S code from Version 1.6 to Version 1.10. 


UPDATE NAME FUNCTION 

BCS17 ALLOWS ALTERNATE INPUT FILE FOR 

AEROELASTIC DATA AND ADDS 
SELECTIVE DEGREE-OF-FREEDOM MATRIX 


CONTROL ADDS 2 SOLUTION TASKS: 

(A) STATIC RIGID ANALYSIS WITH 
CONTROL SURFACE DEFLECTION 

(B) STATIC AEROELASTIC ANALYSIS WITH 
CONTROL SURFACE DEFLECTION 


INVERR ADDS OUTPUT OF AEROELASTIC DATA 

TASKCORR UPDATES DOCUMENTATION AND LOGIC FOR 

CONTROL UPDATE 


INTEQM 


ADD OPTION FOR INTEGRATION OF EQUATIONS 
OF MOTION FOR ZERO AIR FORCES 
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MISCELLANEOUS PROGRAM ENHANCEMENTS (continued) 


UPDATE NAME 
NEWBL 

SHOCKFIX 

I 

i 

VERS110 

WRTWRNS 

* 

WPAFB1 

WPAFB2 


FUNCTION 

REMOVES EFFECT OF SHOCK SWEEP ON 
BOUNDARY LATER CALCULATION 

PROVIDES PROPER SHOCK PLACEMENT 
FOR NEWBL 

UPDATES PROGRAM BANNERS TO VERSION 1.10 
ADDS ADDITIONAL ERROR DIAGNOSTICS 
UPDATES FOR INSTALLATION AT WPAFB. 
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IMPROVED VECTOR I ZATION OF XTRAN3S 


The original version of XTRAN3S was developed using the CDC 7600 computer. 
During the code development period, advanced vector machines, such as the 
VPS-32 and CRAY-1S and later the CRAY X-MP became available. The code was 
adapted to operate on these machines, but only minimal changes to the code 
were performed to take advantage of "implicit" vectorization, i.e., those 
portions of the code that could easily be adapted to vector computation 
through the existing nature of the computational algorithm. 

Later studies with a pilot code version of XTRAN3S, operational on the 
CRAY-1S, showed that at least a factor of two improvements in computational 
efficiency could be achieved by rearrangement of operations to permit a 
larger degree of implicit (or automatic) vectorization. (This represents a 
speed-up factor of almost fi ye compared to an unvectorized or scalar code 
operating on the same machine.) These concepts are generally applicable, 
but not directly transportable, to the VPS-32 version of the code, since 
VPS-32 requires longer vectors than the CRAY-1S to achieve improved 
efficiency when compared with scalar computations. The improved 
vectorization scheme has now been implemented on the CRAY X-MP in Version 
1.10 of XTRAN3S. 

In this section of the presentation the original algorithm will be 
described with respect to its implications for vectorization. Then the 
modifications necessary to achieve a higher degree of vectorization on the 
CRAY X-MP will be discussed. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


ORIGINAL ALGORITHM STORAGE SCHEME (V ersion 1.5) 

In the original version of XTRAN3S, the equations for the X-sweep are 
solved sequentially for each Y mesh point, with a constant Z mesh point 
location. Thus for the default values of Ny = 60, Ny = 20, = 40, 

20 equations of length 60 are solved for each X-Y plane. Then, for the 
same X-Y plane, the Y-sweep equations (60 equations of length 20) are 
solved sequentially in the downstream (increasing X) direction due to the 
presence of the backward spatial or "upwind" difference approximation to 
♦xt* This process is then repeated for the next X-Y plane in the 
increasing Z direction*. The Z-sweep is performed by accessing the data, 
formulating the equations, and solving for each X-Z plane sequentially in 
the increasing downstream direction, with the process then repeated for the 
next plane in the increasing n direction. Since solutions of the Z-sweep 
equations, = (j> n+1 , are dependent on <t> n , and the solution of the Y-sweep 
equations but not on the solutions to the X-sweep equations, <i> need not 
be stored in a three-dimensional array. 
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SOLVE 20 EQUATIONS OF LENGTH 60 
REPEAT 40X 

(EQUATIONS ARE INDEPENDENT) 




Y SWEEP 




SOLVE 60 EQUATIONS OF LENGTH 20 
REPEAT 40X 

(EQUATIONS ARE DEPENDENT-MUST BE 
SOLVED SEQUENTIALLY IN i DIRECTION! 


k 




SOLVE 60 EQUATIONS OF LENGTH 40 
REPEAT 20X 

(EQUATIONS ARE DEPENOENT-MUST BE 
SOLVED SEQUENTIALLY IN i DIRECTION! 


*3-0 STORAGE REQUIRED ■ 3 LEVELS - 144,000 
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MODIFIED ALGORITHM STORAGE SCHEME (Version 1.10) 


For the CRAY X-MP, the scheme illustrated in this figure has been adopted. 
For the X- sweep, data is accessed for each Y-Z plane (rather than each X-Y 
plane), and vectorized solution is performed in the Z direction, i.e., each 
step in the solution is a vector operation of length 40. This is then 
repeated for the next Y-Z plane, and the intermediate result 5 is stored in 
a three-dimensional array. For the Y sweeps, the data is accessed in the 
Y-Z planes, and a vectorized solution of length 40 is again performed in 
the Z direction, with the results $ stored in a three-dimensional array. 
Finally, the Z-sweeps are performed as a vectorized solution of length 20 
in the Y direction, and the advanced values of the potential <f> n+ l are 
stored in place of the values <t> n at the previous time step. It should be- 
noted that four three-dimensional "levels" of storage, 4> n , 4> n_1 , 5> and $> 

are required for this scheme, compared with three for the original scheme. 
In addition, however, nine additional three-dimensional arrays for vector 
equation coefficients and right-hand sides have been stored to improve 
efficiency and decrease the amount of re-calculation required. For the 
default mesh, the amount of storage for three-dimensional arrays has 
Increased from 144,000 to 624,000. The total storage requirement has thus 
Increased from about one-half million words to over one million words. 
Because of the larger amounts of storage available on both the CRAY X-MP 
and VPS-32, compared to the CDC 7600, this has caused no difficulty. 




ft* -w* 

♦ # 


SOLVE 1 VECTOR EQUATION 
(VECTOR LENGTH 40) 
REPEAT 20X 



SOLVE 1 VECTOR EQUATION 
(VECTOR LENGTH 40) 
REPEAT 60X 
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SOLVE 1 VECTOR EQUATION 
(VECTOR LENGTH 20) 
REPEAT 60X 


*3-0 STORAGE REQUIRED - 4 LEVELS + 9 COEFF ICIENTS - 624.000 


ORIGINAL PAGE IS 
OF POOR QUALITY 


i 
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VECTORIZATION OF NASA WING-BODY MODIFICATION 


This figure shows the improvements that can be achieved by careful attention 
to the use of vectorization concepts. The NASA-Langley developed wing-body 
modification update was installed in both the standard Version 1.5 and 
improved Version 1.10. Due to the use of conditional branching within 
computation loops (DO-loops), vectorization was inhibited. Version 1.5 
took more than twice as long to run a wing-body case as a wing-alone case 
with the same dimensions, and Version 1.10 took more than three times as 
long. (The unmodified Version 1.10 runs more than twice as fast as Version 
1.5). When the conditional branches were constrained to the outside of 
computation loops, permitting vectorization to occur, considerable speed-up 
was achieved. A wing-body case can now be run using Version 1.10 with only 
an increase of 4 % over the cost of a wing-alone case. These constraints 
are applied in all the modification of XTRAN3S described below. 


VERSION/CASE WING ONLY WING-BODY 


1.5 

1.5 + W-B 

1.10 

1.10 + W-B 
1.10 + W-B 
(VECTORIZED) 

(NORMALIZED 


1.0 

- 

1.05 

2.18 

.43 

- 

1.69 

1.744 

.43 

.45 


CP SEC/ITERATION) 
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XTRAN3S ALGORITHM IMPROVEMENT 


The XTRAN3S code solves the modified transonic small disturbance equation 
shown in this figure. After transforming the wing planform in Cartesian 
Coordinates (X, Y, Z) to a rectangular region in computation coordinates 
U, n. c), the equation appears as shown below. 

The underlined terms originate with the spanwise (Y-dependent) or cross 
(X-Y dependent) terms of the original equation, Ga>y2 and 0y+H<t> X( w. The 
non-underl ined terms depend only on original streamwise or normal 
derivatives. In the XTRAN3S algorithm, the original spanwise term G$ y 2 was 
split into streamwise and transformed normal components. These terms were 
then treated separately with the streamwise component treated implicitly, 
and the normal component treated explicitly. 


MODIFIED TRANSONIC SMALL DISTURBANCE EQ. 


Uv. 


h < * (e% - r<P« * 




— o 


TRANSFORMING TO 


f r iT^yJ ; j? • y J J ' P j 'Tr - t 


7 f 
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The equation can be rewritten as shown by recollecting terms by 
coefficients of like-terms. Here a and b are the^original (x-direction) 
streamwise terms of the governing equations, and a, 6, etc., contain 
coefficients based on the spanwlse terms of the equation. It has been 
previously noted that a time-accurate, fully conservative, implicit 
treatment can only be performed on the terms 

fe (a « + ) and fj (j ♦„) 

The remaining terms can be differenced in a fully conservative explicit 
fashion, or in a partially conservative implicit fashion. In the original 
XTRAN3S algorithm, only the cross- term (G<j> n 2 ) is so treated; the 
remaining terms are treated explicitly, in the present study, several 
different combinations of implicit and explicit difference treatments have 
been investigated. By examining the coefficients it can be seen that 
certain terms which were treated explicitly in the original XTRAN3S 
algorithm can be treated implicitly without loss of accuracy (i.e., without 
losing the conservation form of the equation). It has also been noted by 
Batina' that only th§ original streamwise terms, represented by the 
coefficients a and b should be included in the type-dependent difference 
scheme or in the definition of the shock point operator. This treatment of 
these terms corresponds directly to the original 2-dimensional method of 
Ballhaus and Goorjian known as LTRAN2. 


RECOLLECTING TERMS: 

$(-&"*) 

t<fs ' Vj * 

' fVt 4 

* ^ ) * ° 

WHERE 

a r Ef x ir * (***») i'jVj £ * 

& - te+iOsy , c * f yAx d = <r 

* ’ Vf, f • HI, J - ±4 W(p 5 4- ± 
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ORIGINAL ADI ALGORITHM 


The transformed streamwise sweep for the original (Version 1.5) XTRAN3S 
algorithm is shown in this figure. 

Here Ds is a type-dependent mixed difference operator based on the sign of 
(a n + 2 $5 is a central difference operator, On is a mixed difference 

operator based upon the sign of 2G s (6 n <j>n) . 

C - SWEEP: 

) - »r[* 

+ + 9<P'J 

+ i*„r + **x’ + s, r ‘ 

WHERE 

x” = t( P,) + H 

4 4 s y? x < 

Y" * VrX + H 4 H Oj" 
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MODIFIED ADI ALGORITHM 


The transformed streamwise sweep for the modified (Version 1.10) XTRAN3S 
algorithm is shown in this figure. Here it can be seen that all but one of 
the explicitly-treated streamwise derivative terms («^X n ) in the original 
algorithm can be treated implicitly. The remaining cross-derivative terms 
6 ^( 4 > n ) or 6 n (<t>£) must still be treated explicitly. 


5 - SWEEP: 

s5 > ( tf. ) • • W,] 

* s,[ rfWj « «<?;<?, ] 

"2. '-J 

4 i 4 + «> y " 

WHERE 

x ” = *»/*•« % 

y" = %,% ♦ * f, <?; + h ( p;tf 
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The transformed spanwise and vertical sweeps are unmodified in the new 
algorithm since these solutions are fully implicit. 


n - SWEEP: 

B 


c - SWEEP: 

A 

5* 


( t ^)- i'ldw-iW) 

Q(Cjf) - s f ?' ) 


) + si c < a_j_ 

At*- J ' / 


+ i-± r„ (<?'”-<?) 


(THESE SWEEPS UNMODIFIED) 
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A-6E TEST PROBLEM 


The modified ADI algorithm has been tested on the wing planform for the 
A-6E attack aircraft with a moderate sweep and taper. The computation time 
step could be increased by a factor of ten with numerical stability 
maintained. For a more severely swept and tapered planform (the F-5 wing) 
an increase of more than fifteen has been reported. 



A-6E Planform In XTRAN3S Calculations 
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This figure shows the resulting steady pressure distributions for the original 
and modified ADI algorithms. In the modified algorithm the shock has moved one 
grid point upstream. Otherwise results are essentially identical. 




WING-PYLON STORE TRANSONIC ANALYSIS 


This section of the presentation will describe the development and current 
status of a pylon/external store or nacelle capability in XTRAN3S. 


Wing - Pylon - Store 
Transonic Analysis 
With XTRAN3S 
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A-6 WING WITH STORES 


Flutter characteristics of military aircraft may be severely affected by 
the presence of external stores. These stores may affect the flow field 
characteristics as well as the structural dynamic characteristics of the 
aircraft. Flutter speeds with carriage of external stores are usually 
lower than for the clean wing and may require severe operational 
restrictions. This figure shows a schematic of the wing of the Navy A-6E 
attack aircraft. As a part of the flutter clearance program for the newly- 
designed graphite-epoxy A-6 Replacement Wing, Boeing performed wind tunnel 
tests on both rigid and aeroelastic transonic models. The A-6 wing is 
shown here with 2 -400 gallon external fuel tanks mounted on pylons under 
the wings at two inboard stations (the third station is used to mount a 
small air-air missile such as the AIM-9L Sidewinder). The wing is 
moderately swept and tapered. Simulations conducted to date with XTRAN3S 
did not include the inboard strake, fuselage, or nacelle. 



ORIGINAL PAGE IS 
OF POOR QUALITY 

















As seen in this figure, aerodynamic interference due to the presence of the 
stores has a significant effect on the sectional lift distribution of the 
wing. (Symbols represent data from a static rigid test). Also shown is a 
prediction of the static lift curve slopes by XTRAN3S for the clean wing. 
Flutter results using the clean wing XTRAN3S aerodynamics and the clean 
wing test data as correction factors in an interfering modified strip 
theory analysis are essentially the same. 


A *6 WING M « 0.87 


•16 \ 


XTRAN3S 


O CLEAN 
□ 3 STORES 



FLUTTER TEST COMPARISON - ONE STORE 


In this figure, analysis of the A-6 wind tunnel flutter model tested in 
NASA's Transonic Dynamics Tunnel for both the "clean wing" and "stores on" 
aerodynamics correction factors (section lift curve slopes and aerodynamic 
centers) are compared with measured flutter dynamic pressures. The 
analysis for both sets of corrections shows the typical "transonic bucket" 
as a function of Mach number. The model test with an actual stores 
representing a single 300 gallon external fuel tank showed essentially no 
effect of Mach number - a completely unexpected result! An effort to 
explain this phenomenon, a "pencil store", or simulator with identical mass 
and inertia properties but minimum aerodynamic interference was also 
tested. Here the results showed the expected Mach effect of decreasing 
flutter speed up to a critical Mach number. A "no-flutter" run (not shown) 
at higher speeds demonstrated the classical transonic bucket behavior. 

Thus the importance of wing-store aerodynamic interference and potential 
aerodynamics of the store itself were demonstrated. 



ORIGINAL PAGE IS 

OF POOR QUALITY 




ORiiGtNAL PAGE- S'b 

OF POOR QUALITY 


FLUTTER TEST COMPARISON - TWO STORES 

This figure shows a similar comparison of analysis and test results for two 
300 gallon tanks and pencil stores. Although the measured flutter boundary 
shows a much higher Mach number for the same dynamic pressure at flutter 
(or vise versa), the bucket behavior for the pencil stores is again 
confirmed. The actual stores show no appreciable effect of Mach number on 
flutter. One principal difference between the actual stores and the 
simulators was the rapidity of the onset of flutter. With the actual 
stores, flutter was approached quite slowly with fairly long times to 
double amplitude (5-10 seconds) allowing good confirmation of actual 
flutter behavior without risk to the model. On the second flutter point 
shown with the simulators (at a q of 155 psf) violent flutter was 
encountered (.1 sec to double amplitude) and the model was destroyed. 



FLUTTER ANALYSIS, 

TWO 300 GALLON TANKS ON PYLONS 4&5 
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(SIDE OF BODY) 


INITIAL XTRAN3S STORE-PYLON SIMULATION 

This figure shows a portion of the grid (Y-Z plane) used in the internal 
XTRAN3S simulation of the A-6 wing-pylon-store geometry. The pylons are 
located in planes midway between spanwise grid lines. The stores are 
simulated as a single line, with stores boundary conditions being satisfied 
at the four neighboring grid points. Shown in this figure are the maximum 
diameters of the 400 and 300 gallon tanks in relation to the wing and an 
approximate side of body location. The rigid model was tested with the 
store configuration shown plus a Sidewinder at the outboard station (which 
has not been modeled). 


INITIAL XTRAN3S STORE-PYLON SIMULATION 
(MING LOWER SURFACE) 


A-6 CLEAN NING 


This figure and the next show the measured pressures on the wing upper and 
lower surface at a location just outboard of the inboard pylon station. 
Shown are the clean wing and pylons only cases. Also shown are the 
pressures predicted by XTRAN3S including surface boundary layers. Other 
than difference at the leading edge and slight pressure differences on the 
lower surface, the comparison is good. 


A-6 CLEAN MING 


M.- .87 - .0 

Tetr XiHAkJiS 

< r 

% 

□ 


O' 


.07* ./ 2.0 



i 


l 
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A-6 PYLON + STORES 


This figure shows the comparison of measured and predicted pressures for the 
pylon-stores case. Although there is additional interference due to this 
"slender-body" modeling of the store, the comparison is neither 
qualitatively nor quantitatively good. Thus is was seen that a more 
complete modeling of the store is required. 


M ' . 87 #.= o ° 


>1 = . 33 


7£st XTfa^ZS? 



INTERFERENCE SHELL GEOMETRY 


This figure shows the effective geometry of an "Interference shell" or 
surface on which the boundary conditions are satisfied, surrounding the 
store. In the manner of Batina's Wing-Fuselage modification, slender-body 
corrections are used to modify the local surface boundary condition applied 
on the shell. 



ORIGINAL PAGE IS 
OF POOR QUALITY 
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(SIDE OF BODY) 


XTRAN3S STORE-PYLON SIMULATION 


i 




This figure shows the grid geometry for the store-pylon simulation with 
interference shells. 




A-6 PYLONS + STORES (INTERFERENCE SHELL) 


In this figure the XTRAN3S prediction with the interference shell is 
compared with measured data. The analysis predicts a much stronger shock, 
located closer to the trailing edge of the wing, than the data indicated. 
(The boundary layer calculation for this case was unstable; thus the 
result shown is an inviscid calculation.) 

(At the symposium, the possibility of flow field separation causing the 
behavior shown was suggested by Joe Giesing of McOonnel 1 -Douglas. This has 
been simulated and the results are shown in the next figure.) 



42 


o iJ 



A-6 PYLON + STORES WITH AFT SEPARATION 


Flow separation on the aft portion of the stores was simulated by use of a 
constant store radius and zero slope for the portion of the store aft of 
the maximum radius. This geometry modification has considerably improved 
the correlation with measured data, as shown in this figure. This 
calculation is shown with the boundary layer present. (This figure was not 
shown at the symposium.) 

I 

i 



OM'GiNAL PAGE-! IS 
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SUMMARY 


In summary, modifications to the XTRAN3S code, implemented in Version 1.10, 
have provided a speed-up of twenty to thirty times over the previous 
released Version 1.5. 

The pylon-store modification, under development, has shown that the 
presence of pylons and stores can be simulated although attention to 
geometric and flow field details is required. 


• IMPROVED VECtORIZATION and ALGORITHM 
MODIFICATIONS PROVIDE 20-30X SPEED-UP 
IN XTRAN3S VERSION 1.10 


• PYLON-STORE MODIFICATION DEVELOPED 

• SIMPLE SLENDER BODY INSUFFICIENT 

• INTERFERENCE SHELL OVERPREDICTS STORE 
INTERFERENCE ON WING 



FUTURE EFFORTS 


Future developments of the XTRAN3S code, to be initiated shortly under 
AFWAL sponsorship, are shown In this figure. Some comparisons with 
NASA- Lang ley's CAP-TSD code will be made at a later date. 


• CONTINUE DEVELOPMENT ON PYLON- STORE 

MODIFICATION 

• INCORPORATE WING-BODY AND PYLON-STORE 

MODIFICATIONS IN SINGLE VERSION 

• ADD SUPERSONIC AND FOURIER ANALYSIS CAPABILITIES 

• CONTINUE EVALUATIONS WITH TRANSPORT AND ATTACK 

AIRCRAFT AEROELASTIC MODELS 



